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Consider a spherically symmetric spacetime generated by a self-gravitating massless scalar field <j> 
and let ip be a test (nonspherical) massless scalar field propagating on this dynamical background. 
Gundlach, Price, and Pullin [3| computed numerically the late-time tails for different multipoles 
of the field ip and suggested that solutions with compactly supported initial data decay in accord 
with Price's law as t~^ 2l+ ^ at timelike infinity. We show that in the case of the time-dependent 
background Price's law holds only for I — while for each I > 1 the tail decays as t~ f2e+2 \ 
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The Einstein-massless scalar field system 

Ov 

G a0 = 8tt ( V a 4>Vfi4> - ^(V^V'V) , g a N a V^ = 0, (1) 

<N : V 2 y 

O \ restricted to spherical symmetry, has been serving as an important theoretical laboratory for the investigation of 
nonlinear gravitational phenomena in a rather simple 1 + 1 dimensional setting. For this system Christodoulou proved 
| that a generic spherically symmetric solution settles down asymptotically either to Minkowski spacetime (for small 
[ — . data) [l| , or to a Schwarzschild black hole (for large data) [2j • The first reliable numerical simulations of the late-time 
t— I ' asymptotics of this relaxation process have been done by Gundlach, Price, and Pullin (GPP) 3]. They found that, 
regardless of the endstate of evolution, the scalar field develops a tail which falls off as t~ 3 near timelike infinity (for 
' q ' compactly supported initial data). 

q-i In a recent paper Q we revisited this problem to emphasize that the asymptotic convergence to a static equilibrium 
(Minkowski or Schwarzschild) is an essentially nonlinear phenomenon which cannot, despite many assertions to the 
contrary in the literature, be properly described by the theory of linearized perturbations on a fixed static asymp- 
totically flat background (Price's tails d, Q). This is particularly evident for dispersive solutions which asymptote 
Minkowski spacetime. In that case the quantitative characteristics of the tail (the decay rate and the amplitude) can 
be obtained using nonlinear perturbation expansion [Jj. Since some details of this formal calculation will be needed 
' below, let us now briefly summarize it. In the parametrization 

i> ' _, 

ds 2 = (l- (_ & mt,r) d1 2 + dr 2^ + r 2 {M 2 + sm 2^ 2) ^ (2) 

(N 

the system ([1]) takes a particularly convenient form (below an overdot denotes d/dt and a prime denotes d/dr) 

ON ' 

ml = 2nr(r -2m){<j)' 2 +e~ 2/3 2 ) , (3) 

m = An r(r — 2m)(j) 4> , (4) 

& - ^r~^ (5) 



r(r — 2m) 

= ^(r 2 eV)'- (6) 

Consider small and compactly supported initial data (<f), <f>)t=o = ( £ f( r ), e 9{ r ))- Then, up to the order 0(e 3 ), we have 

4> = £0i + e 3 3 , m = e 2 m 2 , (5 = e 2 f3 2 , (7) 
where <f>\ satisfies the flat-space radial wave equation 

□fc := 01 - 0? - Vx = 0, (0 1 ,0 1 ) f=o = (/(r),.g(r)), (8) 
r 

while the second-order perturbations of the metric functions satisfy 

r4 = 2 7 rr 2 (0 2 +0f), 0' 2 = ^ . (9) 
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Solving equation ([5]) and then © we get for t > R (where R is the radius of support of initial data) 

a(t — r) , 

Mt,r) = (10) 

m 2 (t,r) = 2tt I 2 / a /2 (s)ds- gJL_LZ I ; (H) 
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4tt j 
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.;,(/./•) = | - / a'\s)ds + 2 j ^ds- j -0^ds\ . ,12) 

t— r i — r 

where the initial-data-generating function a(u) vanishes for |u| > R. The third-order perturbation of the scalar field 
4>3 satisfies the inhomogeneous wave equation (with zero initial data) 

= 2/? 2 ^i + Mi + PWi =■ S(t, r) . (13) 
The source S(t,r) is already known from (10-12) so we can use the Duhamel formula 

t t+r-T 

Mt,r) = ^Jdr J P S(T,p)dp, (14) 
\t-r-r\ 

to obtain the asymptotic behavior for large retarded times 

<f>(t,r)~e 3 Mt,r)~ (t 2^ )2 ■ ^ = -2 5 tt J a(u) J {a' {s)) 2 ds du . (15) 

— OO u 

We refer the reader to Q for more details about this calculation and numerical evidence. 

After this introduction, we are ready to discuss an interesting model for investigating linear nonspherical tails on 
a fixed dynamical background. This model, proposed by GPP 3], involves a nonspherical test massless scalar field 
ip which propagates on the spacetime @ generated by the self-gravitating field (j). Since the dynamics of ijj is linear 
and the background is spherically symmetric, one may decompose tp into spherical harmonics 

il>(t,r,0,<p)= Yl i>em(t,r)Yr(#,<p), (16) 

£>0,\m\<£ 

and analyze the evolution of each multipole separately 

GPP conjectured 1 that for compactly supported initial data the multipoles have the tail i/j£ m (t,r) ~ i _ ( 2£ + 3 ) at 
timelike infinity (t — > oo at a fixed r), in accord with Price's law on a fixed Schwarzschild background, even though 
the actual spherical background is time dependent and its Bondi mass decreases (to a positive value in the collapsing 
case or to zero in the dispersive case). It seems that GPP's conjecture was based more on belief than numerical 
evidence, because for the first few multipoles the following power-law exponents of the tail were reported numerically 
(see Fig.12 in [|): -2.77 (I = 0), -3.95 (£ = 1), -5.94 [t = 2), and -8.34 ({, = 3). 

The purpose of this note is to point out that in this model (and for other time-dependent backgrounds) Price's law 
(i.e., £~( 2£ + 3 ) decay) holds only for 1 = 0, while for I > 1 the power-law exponent of the tail is equal to — (2£+ 2) (as 
was clearly indicated by GPP's own numerics). To show that we shall compute the late-time asymptotic behavior of 
ipZmitjr) (for smooth initial data compactly supported in a ball of radius R') along similar lines as described above 
for the field <fi. The perturbation expansion has the form ip£ m = ipo + e 2 ip2 + ■ • ■ , where for convenience of notation 
we dropped the multipole indices on iterates. At the zero order we have 

2 £(l + l) 



1 GPP considered the characteristic initial value problem while we are studying the Cauchy problem, however this difference do not affect 
the asymptotics of tails. 
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which for t > R' is solved by 



1 JL, (2l-k)\ bW(t-r) 



r ^ k\{£-k)\ (2ry- k ' 



(19) 



where the initial-data-generating function b(u) vanishes for |u| > R' (the superscript in round brackets denotes the 
fc-th derivative). At the second order we get 



□ W V>2 = 2/3 2 ^o + Mo + PWo ~ 2£(£ + 1)m2 ^ =: Se (t, r) , ty> 2 , ^)t=o = (0, 0) . 



(20) 



Substituting (11), (12), and p^|) into the Duhamel formula (where Pe(x) is the Legendre polynomial of degree £) 

t t+r-T 



|t-r- t| 
we get (for large retarded times) for i = 0: 



r 2 + p 2 - (t - t)' 
2rp 



Se(r, p)dp, 



(21) 



V»oo(t,r) ^ £ 2 ^2(t,r) 



e 2 B t 



+oo +c>o 

B = -2\ { b(u) { {a'{s)) 2 dsdu 



-OO It 



(t 2 - r 2 ) 2 ' 

and for I > 1 (see 0] for technical details of calculation in the I > 1 case): 

£ 2 £^ „ , _,2 £+3 ^!7r + Ff£ 2 (£- 1) 



(22) 



tpim(t,r) ~ e ip 2 (t,r) 



B^i-^^TVT- I l^--^a 2 (u)b"(u)~2£ 2 (a'(u)) 2 b(u)) du. 



(t2-r 2 Y+ 1 ' ~* v *' 2^+1 J V 2£-l 

— OO 

The numerical verification of these formulae is summarized in Figure 1 and Table 1 
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FIG. 1: Left panel: The difference between the local power index and the theoretical prediction: n(r = l,t) — 3 for I = and 
n(r = l,t) - 2(£ + 2) for I = 1, 2, 3 , as a function of 1/t. Right panel: The log-log plots of \tp{t, r = 1)| for £ = 0, 1, 2, 3. Both 
panels correspond to initial data generated by a(x) = exp(— x 2 )/^/2tt (for 4>i) an d b(x) = x 2 exp(— x 2 ) (for ipo), with e = 2~ 8 . 



We fit our numerical data with the formula 

ip(t, r) = Bt" 1 cxp (c/t + d/t 2 ) , (24) 

which gives the local power index (LPI) [8] 

n{t, r) := -tip(t, r) /ip(t, r) = 7 + c/t + 2d/t 2 . (25) 

Our fitting procedure proceeds in two steps. First, from the local power index data on the interval 1/t < 1/50 (the 
left panel of Figure 1) we fit 7, c and d in (|25[) . Next, having determined 7, c and d in this way, we fit B in (|24[) from 
tp data on the interval 50 < t (the right panel of Figure 1). The results of this procedure are given in Table 1. 
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Numerics: LPI data 


Theory (third order) 


Numerics: ip data 


c | d | 7 


7 | B 


B 


£ = 




-3.068 


6.056 


3.000 


3 


-5.296 x 10~ 4 


-5.293 x 10~ Y 


2 -iu 


-3.064 


5.877 


3.000 


3 


-8.474 x 10~ 4 


-8.468 x 10~ b 




-3.064 


5.867 


3.000 


3 


-1.356 x 10~ 4 


-1.355 x 10~ 4 


1= 1 


2 -i2 


16.73 


-134.7 


4.008 


4 


1.084 x 10 - '' 


1.145 x 10"' 


2 -iu 


16.73 


-134.7 


4.008 


4 


1.735 x 10 -0 


1.833 x 10" b 


2-8 


16.73 


-134.7 


4.008 


4 


2.776 x 10~ b 


2.932 x 10~ b 


£ = 2 




15.87 


-139.4 


6.005 


6 


-6.940 x 10"' 


-7.193 x 10~ Y 


2 -iu 


15.87 


-139.4 


6.005 


6 


-1.110 x 10~ 5 


-1.151 x 10~ b 




15.87 


-139.4 


6.005 


6 


-1.777 x 10" 4 


-1.841 x 10~ 4 


£ = 3 


2-1^ 


13.28 


-110.7 


8.012 


8 


6.023 x 10"° 


6.525 x 10" b 


2 -iu 


13.29 


-111.0 


8.012 


8 


9.637 x 10~ b 


1.044 x 10~ 4 


2- 8 


13.31 


-111.4 


8.012 


8 


1.542 x 10~ a 


1.669 x 10 -3 



TABLE I: The comparison of analytic and numerical decay rates and amplitudes of the tails at timelike infinity for I = 0, 1, 2, 3, 
for initial data generated by a(x) = exp(—x 2 )/\/2iv (for cf>i) and b(x) = x 2 exp(—x 2 ) (for ^o), at r — 1. The theoretical 
prediction is B — e 2 Be, with Bo given in (|22p and Bi given in (123[l for t > 0. Fits were made on the interval 50 < t < 200. 

It is instructive to compare the tail (|23| with the tail of a massless scalar field propagating on a fixed asymptotically 
flat static background. The latter can be readily obtained from the Duhamel formula (|2"Tj) applied to the source (|20|) 
with m2 = M and /3a = — 2M/r, where M is the total mass. The result (valid for all i) reads 

il> tm (t,r) - {t2 C _ e l, 2 y +2 , Q = -M2 i+3 (£+l)\ I b(u)du. (26) 

— oo 

This is the celebrated Price tail 5] (as far as we know, first obtained in the form (|26[) by Poisson 0). It is worth 
stressing that the formula ([26]) yields a good approximation for the amplitude of the tail provided that both an observer 
and initial data lie in a weak field region where M/r is small. Note that for I — the formula (|22|) takes the form 

(|26)> when the total mass M is replaced by the weighted average over the Bondi mass M(u) = 4ir J (a'(s)) 2 ds. For 

u 

I > 1 Price's tail decays by one power faster than that in (|23| which, on a technical level, is due an extra cancelation 
in the integration by parts of Duhamcl's formula. 
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